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20.  ABSTRACT  (firtAii  «n  r*  rm*mm  if*  it  ntfMiy  mod  idmrui  tj  br  block  m mtlmt) 

The  self-diffusion  coefficients  of  ^'C  in  NbC..  sinqle  crystals  have  been 
measured  as  a  function  of  composition  in  the  temperature  range  1900-231 SK, 
and  can  be  represented  by  the  expressions 
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The  lower  values  of  the  activation  energy  and  the  pre-exponential  term  in 
NbC  „  are  attributed  to  a  change  in  the  path  of  C  mass  transport  from  that 
of  an  octahedra  1  - totrahedra  1  -octahedra  1  mechanism  in  Uhl ,  ; .  H,,lt  and  *'hCo.«?/,  to 
one  involving  a  r-mef.u!  di vacancy  mechanism.  3  ho  eltect  of  lattice  geometry 
and  the  electronic  charge  distribution  on  the  diffusion  mechanism  is  also 
discussed. 

The  self-diffusion  of  ',lM>  in  single  cyrstals  of  NbC0>  ace)  NbCo.rm  and 
Nbro  76(,  has  also  been  studied  in  the  range  of  23/0-2660K.  Tfio  diffusion 
coefficients  are  composition  independent  and  can  be  described  by  the  expression : 

D^v,  =  ( 4 . 54  *  j  • )  exp (  -( 1 4 0 . 0 .•  2 . 4Kca  1  /mo  1  e )  / RT ) cm  /  s ec 

Vi  analysis  of  the  results  indicates  that  Nh  diffuses  by  an  (0-0)  mecnanisni, 
.just  as  a  pure  metal  diffuses  in  a  f.c.c.  lattice,  wherein  the  atom  migrates 
from  its  lattice  position  directly  to  an  analogous  vacant  site.  As  this 
process  involves  the  energies  of  both  migration  and  Nb  vacancy  formation, 
the  slower  diffusion  rates  of  this  species,  relative  to  that  of  C  whose 
mechanisms  of  transport  involve  only  the  energy  of  migration,  are  thus 
explained. 


The  conceptual  design  and  construction  of  an  apparatus  developed  for  de¬ 
forming  materials  by  subjecting  them  to  constant  compressive  stresses  at 
elevated  temperatures  is  alxo  described.  The  results  of  actual  creep 
studies  have  shown  that  all  compositions  undergo  rapid  and  extensive  primary 
creep  even  at  low  stresses  (5-10x10  ’ps  i )  arid  temperatures  ( 1400- 1650 "'C)  at 
the  onset  of  loading.  Steady  state  conditions  occur  only  after  many  hours 


^  r-V  C*  V 
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(days)  have  elapsed.  The  steady  state  exponent  is  between  4.2  and  4.5 Ve 
which  endicates  a  diffusion  controlled  mechanism  of  cre^p. 


DIFFUSION  AND  CREEP  IN  NIOBIUM  CARBIDE  AS  A  FUNCTION 
OF  TEMPERATURE  AID  CARBON  CONTENT 


I  .X^atement  of  Problem 

The  transition  metal  carbides  possess  a  combination  of  some  of  the 
highest  melting  points,  hardness  values  and  mechanical  strengths  of  any 
known  materials.  The  contemporary  interest  in  these  materials  in  structural 
and  cladding  applications  at  temperatures  wher>'  most  materials  would  be 
molten  or  wou.d  seriously  degrade  made  it  imperative  that  the  high  temper¬ 
ature  diffusion  and  deformation  processes  and  their  relationship  to  each 
other  be  understood.  Thus,  the  principal  thrusts  of  this  ARG-sponsored 
program  have  been  the  conduct  of  self-diffusion  studies  of  1  *C  and  "Nb 
in  single  crystals  of  N bC^  9 NbC0^11(  and  NbCo  and  constant 

stress  compressive  creep  research  on  the  same  materials.  Extensive 
characterization  of  these  materials  has  also  been  performed. 

1 1 .  Summary  of  Important  Results 

A .  D if f u s i o n  St u d  i e s 
1 .  Carbon 

Typical  penetration  profiles  for  1  C  diffusion  in  the  single 
crystal  specimens  of  NbCx(x- 0.868,  0.B34  and  0.766)  showed  high  activity 
deviations  from  the  straight  line  extrapolated  from  the  deeper  section  of 
the  profiles  toward  the  'near-surface”  (x'O)  portion  except  those  in  which 
the  trace  was  implanted  by  diffusion  into  the  crystals  at  low  (=007310 
temperatures,  the  residual  tracer  removed  mechanically  and  a  NbCx  disk 
placed  atop  the  crystal  during  the  higher  temperature  diffusion  anneals. 

The  "near-surface"  portion  of  the  profiles  was  investigated  in  this  program 
of  research  and  found  to  be  caused  by  factors  extrinsic  to  the  crystal  itself, 
principally  the  d  layed  entry  into  the  crystal  of  that  portion  of  the  tracer 
which  is  not  initially  in  contact  with  the  crystal.  A  complete  discussion 
of  this  phenomenon  from  the  standpoint  of  this  research  and  that  of  others 
is  presented  in  Ref.  1  of  the  bib! iography .  No  evidence  of  high  diffusivity 
path  contributions  to  the  volume  diffusion  was  noted  in  Hie  profile. 


The  linearity  of  the  relative  activity  vs  x'  ulo+s  in  the  deeper 
sections  of  the  profiles  or  the  total  profiles  derived  from  the  1673  K 
implanted  and  covered  samples  indicates  that  these  sections  satisfy  the 
solution  to  Fick's  second  law  for  concentration  independent  diffusion 

from  a  thin-film  source.  As  such  the  calculated  values  of  the  diffusion 

★ 

coeff icier, t,  D.,  represent  volume  self-diffusion  alone.  The  resulting 

^  A 

curves  produced  from  the  plot  of  the  logarithm  of  these  D^.  values  vs  1/T 
for  each  composition.  Figs,  la-c,  reveal  an  excellent  fit  to  the  classical 
Arrhenius  equation  and,  with  the  use  of  computer-aided  least-mean-square 
curve  fitting,  may  be  expressed  as 

Dc^bCo.8C8^niO])-'(2.59!;;®E)xexp(-100.4'2.2kcal/mol)/RT)cm2/s  (1) 

D*(NbC  )([100])=(7.44.!:^)xexp(-105.0i3.3kcal/mol)/RT)cm?/s  (2) 

C  o ,  a  3  c 

D*(NbC  )([110])=(2.22ti - oC)xlO"  xexp(-76. O' 2 . 7kca l/mol  J/RTcmVs  (3) 

Q  0  .  7  f>  f>  • 

The  significant  difference  between  the  activation  energies  and  pre¬ 
exponential  terms  of  NbC0  766  and  the  other  two  compositions  indicates  that 
at  least  two  different  C  diffusion  mechanisms  are  active  in  NbCx  sinqle 
crystals  in  this  ranye  of  C/Nb  ratios  and  that  the  dominance  of  a  particular 
mechanism  is  dependent  on  the  composition.  Furthermore,  these  results  indicate 
that  the  dominant  C  diffusion  mechanism  in  NbC0  r>(<0  and  NbC0.„?u  is  the  same 
but  different  from  that  operative  in  NbCu  766  .  The  slightly  higher  activation 
energy  in  NbCn  with  increasing  C  deficiency  relative  to  that  of  NbC0 

is  consistent  with  the  higher  melting  point  of  former  material  and  is  also 
found  in  TiCx  (2,  3,  4)  and  ZrCx  (5).  However,  NbC0  7f>(,  has  approximately 
tie  same  melting  temperature  as  NbCQ  and  the  largest  C  deficiency  of 
these  three  compositions  but  a  much  lower  activation  energy  for  diffusion; 
thus,  this  is  further  evidence  for  the  existence  of  two  different  C  diffusion 
mechanisms.  A  discussion  of  the  possible  paths  of  transport  are  presented  in 
tie  following  paragraphs. 

From  a  consideration  of  the  rigid  sphere  NaCl-type  structural  model  in 
which  NbCx  forms  and  the  activation  energies  determined  in  this  research,  there 
ore  three  atomic  paths  for  C  diffusion  which  the  present  authors  consider  to 
be  most  probable.  These  are:  (1)  the  (0-0)  mechanism  wherein  the  f,  atom 
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rit).  ].  Carbon-14  sel f-di ffusion  coefficients  for  the  deeper  sections  of  the 
profiles  plotted  as  a  function  of  1/T  for  (a)  NbC0.Br,R’  (&)  NbCc>B'4, 
and  (c)  NbC0  .;r,( . 
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jumps  along  the  [110]  direction  directly  to  the  nearest  neighboring  vacant 
C  octahedral  site  with  a  jump  distance  &0/J[2,2)  the  (0-T-0)  mechanism 
wherein  the  C  atom  jumps  immediately  to  the  nearest  vacant  C  octahedral 
site  via  an  unoccupied  tetrahedral  site  formed  by  the  Nb  sublattice  with 
an  effective  jump  distance  aQA/2  and  3)  the  C-Nb  transient  "di vacancy" 
mechanism  wherein  the  C  atom  jumps  to  the  analogous  vacant  C  octahedral  site 
by  way  of  a  normally  occupied  metal  vacancy  which  has  momentarily  become 
associated  with  the  C  vacancy.  This  last  mechanism  has  been  discussed 
previously  by  Zagriazkin  (6)  and  Andriyevskii  and  Gurov  (7). 

If  only  the  strain  energy  produced  in  making  the  jump  is  considered, 
the  ease  with  which  the  jump  of  each  diffusion  mechanism  noted  above  can 
be  made  depends  on  the  gap  or  separation  between  the  metal  atoms  and  the 
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relative  sizes  of  the  C  and  metal  atoms.  Using  1 ,54A  as  the  atomic  diameter 
of  C  and  assuming  the  atomic  diameters  of  the  various  inetal  atoms  to  be 
those  of  a  hypothetical  fee  metal  lattice,  the  principal  investigator  has 
shown  that  the  C  and  metal  atoms  are  in  contact  along  the  [100]  direction; 
however,  gaps  occur  between  the  metal  atoms  in  the  [110]  directions  having 
a  size,  ?no,  for  NbC  of  0. 231  A-  Since  this  size  is  much  smaller  than  the 
size  of  the  C  atom,  a  large  lattice  strain  or  acconniodating  movement  of  the 
two  Nb  atoms  is  required  in  order  for  the  former  species  to  jump  directly 
from  one  sublattice  site  to  an  adjacent  vacant  C  site  and  therefore,  by 
inference,  for  the  (0-0)  mechanism  to  be  operative. 

By  comparison  the  shorter  [111]  jump  which  would  be  employed  in  the 
(0-T-0)  mechanism  would  have  the  C  atom  passing  directly  through  a  much 
larger  triangular  gap  in  the  (111)  plane,  gin  -  ar,d  into  a  tetrahedral  hole 
at  the  (‘4,  4,  4)  position  in  the  carbide  lattice.  For  NbC,  gm=0.72CA. 

In  essence,  two  jumps  are  involved  in  going  to-and-from  the  tetrahedral 
site;  however,  from  a  strain  energy  of  migration  viewpoint,  of  the  two 
mechanisms,  the  latter  would  have  the  lower  energy.  The  C  vacancy  formation 
component  of  Q  would  be  the  same  in  both  cases. 

However,  from  a  consideration  of  the  relative  volumes  of  the  C  and  Nb 
atoms,  it  is  obvious  that  the  energy  of  motion  component  in  the  activation 
energy  in  the  C-metal  transient  "divacancy"  diffusion  mechanism  will  be 
much  less  than  for  the  motion  of  the  C  atoms  in  the  above  two  mechanisms; 
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since,  in  the  latter  processes,  the  C  atom  has  to  forcefully  miqrate 
between  Nb  atoms.  Therefore,  it  is  concluded  from  a  consideration  of  the 
rigid  sphere  model  and  a  comparison  among  these  mechanisms  that  the  energies 
of  C  migration  in  the  various  mechanisms  decrease  in  the  order  (0-0), 

(0-T-0) ,  C-Nb  di vacancy.  As  the  Q  of  NbC„.7fc  is  less  than  that  of  the 
other  two  compositions,  the  logical  assignments  from  these  points  of  view 
are  the  divacancy  path  to  N b C 0 _  ,(ifi  arid  the  (0-1-0)  path  to  the  other  two 
compositions.  However,  this  does  not  explain  why  the  transport  mechanism 
should  change  at  the  lower  C  concentration  or  why  all  compositions  do  not 
follow  the  same  lowest  energy  path.  In  order  to  understand  this  phenomenon, 
one  must  also  consider  the  changes  in  electronic  charge  densities  in  the 
vicinity  of  the  proposed  avenues  of  transport. 

In  the  hard  sphere  model  employed  above,  all  vacancies  have  the  same 

properties  regardless  of  their  concentration.  However  Williams  [8]  has 

reported  that  the  scattering  of  the  electrons  by  vacancies  in  the  non- 

stoi chiometric  crystals  of  Ti Cx(x-0. 8-0. 97 )  was  not  linearly  proportional 

to  the  vacancy  concentration  and  thus  implied  that  the  C  vacancies  do  not 

have  the  same  properties  in  the  different  compositions.  The  same  is 

believed  to  also  be  true  for  NbCx-  More  recently.  Chad i  and  Cohen  [9], 

in  their  theoretical  research  on  the  charge  densities  in  NbC,  have  shown 

that  in  going  from  NbCjto  NbC„  .  1)  the  number  of  valence  electrons  is 

reduced  from  9  to  8  and  2)  the  magnitude  of  the  charue  in  the  peak  of  the 

d  component  of  the  d  states  becomes  nearly  zero  in  the  latter  material, 
xy  f 

Because  of  the  extension  of  the  d  eleciron  orbits  from  the  neighboring 

xy 

Nb  atoms  into  the  metal  vacancy,  the  C  atoms  diffusing  through  the  metal 

vacancy  would  have  a  barrier  unless  the  electron  density  in  the  d  peak 

xy 

were  close  to  zero.  Thus,  the  feasibility  cf  the  C-Nb  transient  divacancy 
mechanism  is  enhanced  as  the  C/'Nb  ratio  decreases  to  a  value  close  to  0.75. 
Such  evidence  lends  considerable  support  to  the  choice  of  the  divacancy. 
mechanism  for  C  in  NbC,  ....  as  noted  above,  but  the  theoretical  treatment 
also  provides  information  as  to  the  choice  of  a  path  of  transport  for  the 
other  two  compositions. 

The  presence  of  the  large  number  of  structural  C  vacancies  in  all 
compositions  would  indicate  that  the  activation  energy  for  the  (0-T-0)  and  the 
transient  divacancy  mechanisms  would  consist  essentially  of  the  C  migration 
energy  in  the  former  and  the  energy  of  Nb  vacancy  formation  in  the  latter 
path  of  transport. 
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From  vapor  pressure  data,  Storms  et.al.  [10]  have  calculated  the  Nb 
vacancy  formation  energy  for  NbC0.75oto  be  60.49  kcal/mol.  Given  the 
fact  that  there  is  no  doubt  some  small  energy  of  C  migration  associated 
with  the  di vacancy  mechanism,  this  value  for  the  formation  of  the  Nb 
vacancy  is  in  line  with  76.0  kcal/mol  Q  value  found  in  NbC0  766  and 
supports  the  conclusion  that  this  mechanism  is  operative  in  this  material. 

2.  Niobium 

The  95Nb  diffusion  coefficients  calculated  from  the  individual 
profiles  using  tne  thin  film  solution  to  Fick's  second  law  obey  closely 
the  classical  Arrhenius  equation  for  all  three  compositions  of  NbCx.  With 
the  use  of  computer  assisted  least-mean-square  curve  fitting,  these 
coefficients  may  be  expressed*  as: 

D*b  (NbC0.e.8)([nO])  =  (8.99t1i;^)exp(-(143.5±3.8fCcal/mole)/RT)cm2/sec  (5) 

DNb(NbC0.83J([lQ0]H2.00lJ- ^)exp(-(135. 7±3.2Kcal/mole)/RT)cmz/sec  (6) 

°Nb(NbCo  766)([110])+(6.6lt1S; 7i|)exp(-( 142. 0±6. 3Kcal/mole)/RT)cm2/sec  (7) 

The  95Nb  self-diffusion  coefficients  in  these  three  compositions  are 

nearly  equal  to  one  another  and  thus  appear  to  be  essentially  independent  of 

the  C/Nb  ratio.  These  results  are  similar  to  those  reported  by  Sarian  (11) 

for  ** **Ti  in  TiC  .  This  is  also  graphically  illustrated  in  the  combined  plot 
*  x 

of  as  a  function  of  1/T  shown  in  Fig.  2  for  the  three  NbCx  compositions. 
This  resulting  curve  can  be  expressed  as: 

DNb(NbCx)=(4-54-i:7^exp^‘(140-0±2-4Kcal/mo1e)/RT)cmZ/sec  (8) 


*The  notation  at  the  right  of  D.,.  (NbC  )  indicates  the  crystallographic 
direction  of  95Nb  transport;  tne  limits  on  the  preexponential  term  (D  ) 
and  the  activation  energy  (Q)  are  the  probable  errors.  0 
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Kig  2.  N'ioHum-95  seif-<ii(Fusion  cncfticicn’s  vs  the  rt\  »pr<>v.ii  vf 
the  absolute  temperature  for  NK\  <\  -  I1.8M.  0.8 1-;  .iml  0  7f.M 


The  possible  mechanisms  of  migration  of  the  lib  atoms  to  the  nearest  Mb 
lattice  vacancy  in  the  NaCI  structure  of  those  three  compositions  are  the 
same  as  that  described  above  for  C  diffusion.  I"  mu  ter  iu  Is  such  as  Nb(.^ 
wherein  the  metal  atom  contains  a  partially  filled  d  shell  and  the  consti tutents 
are  octahedral  l.y  coordinated,  a  portion  of  the  directions  of  the  d  and  d 
electron  states  in  the  octahedron  around  C  may  be  schematically  portrayed 
as  shown  in  Fig.  3.  Because  of  the  extension  of  the  Nb  d.  electron  states 
into  the  C  vacancy  and  the  resultant  repulsive  forces  which  would  be  felt  by 
a  Nb  atom  in  this  position  and,  more  importantly,  the  much  smaller  size  of  the 
C  vacancy  compared  to  that  of  the  Nb  atom,  the  migration  energy  for  the  latter 
atom  through  the  C  vacancy  will  be  extremely  high.  Furthermore,  the  efficacy 
of  this  divacancy  mechanism  would  be  dependent  on  the  concentration  of  C 
vacancies  which  is  contrary  to  composition-independent,  coefficients  exper¬ 
imentally  determined  for  the  ,5Nb  diffusion.  Thus,  it  is  concluded  that 
the  Nb  atom  does  not  diffuse  via  the  C  vacancy. 

*The  d  states  are  d„  .  dand  d,  and  the  d  states  arc  d  ■  ? 

*  xy  y l  six  i  x  y 


LtJ»bon  t'orrt 


Q  n.,. 


I  '  N»  b'Jf*  vacancy 

y 

(  ig  I  Some  of  lhe  dir.Unmv  of  ihe  d.  .in.!  .1,  'tulev m  Nb.  , 
oU.iheJr.rn  1  he  J,  vi.it, v  of  lhe  Nh  ,i»  me  JueUeJ  toward  the 
C  atom,  and  inti'UU  » ith  Ihe  'r  o.ile.  of  .  .  the  J,  -late,  are 
directed  toward  olher  Nh  atom-  and  I  aim  Nb-Nn  bond- 


From  d  general  consideration  of  the  (0-0)  and  the  (0-T-0)  mechani sms  in 
NbCx,  the  nsNb  migration  energy  depends  on  the  size  of  the  gap  relevant  in 
each  mechanism  whicii  the  Nb  atoms  must  pass  through  to  reach  the  analouous 

o 

vacancy.  Using  the  values  of  1.54  and  2.92  A  for  the  atomic  diameters  for 
C  and  Nb,  respecti vely,  the  gap  size,  gNb[Mij,  encountered  by  the  ’’■‘Nb 
in  diffusing  along  the  (0-0)  path  is  calculated  to  be  1.6  A, which  is 
larger  than  the  0.°4  A  diameter  of  the  tetrahedral  site;  thus,  it  is  con¬ 
cluded  that  Nb  atoms  diffuse  by  the  (0-0)  mechanism  in  NbCy . 

lhe  proper  equation  for  the  Nb  self-diffusion  coefficients  in  NbCx, 
assuming  a  (0-0)  diffusion  mechanism,  can  be  expressed  as 


°Nb=  fN  b  VNbao  exp  ( ( 4Nb+SvNb )  /R  )exp ( -  ( H, „rjl/ Hu|uh)/RT ) 


ir.Nh  vNb' 


where  f^b  is  the  correlation  factor,  v^b  the  Nh  jump  frequency,  a(i  the  lattice 
parameter,  S  and  H  the  entropy  and  enthalpy,  re»p-'C.tively,  for  Nb  atom 
migration  and  5v^b  and  H  the  entropy  and  enthalpy,  respectively,  for  Nb 
vacancy  formation. 

The  composition  dependence  of  the  enthalpy  of  formation  of  the  Nb 
vacancy  in  NbCx  has  not  been  reported.  However,  Kaufman  and  Clouqherty  [12] 


and  Chang  [13]  ^;e  employed  the  Schottky-Wagnor  [14]  model  of  defective 
Hnases  to  derive  the  expression 

A.ft  =  -r-+  +RTln[(l-x)/4]-2RT  In  u  (10) 

M  M 

where  AG  is  the  partial  free  energy  of  formation  for  the  metal  species, 

+  M 

Gc  is  the  free  energy  to  form  the  metal  vacancy,  x  is  the  carbon /metal  atom 
ratio  and  u  is  the  ratio  of  the  vacant  lattice  sites  to  the  total  lattice 
sites  at  stoichiometry  (x=1.0).  The  vacancy  concentration  at  . „oichiometry 
was  assumed  to  be  less  than  1  .  Kaufman  and  Sarney  [15]  further  assumed 

o  o 

that  RT  In  a  =  AH^(0K,x=l .0)  where  is  the  standard  enthalpy  of  formation 
of  the  Nb  vacancy.  Storms  et  al.  [16]  have  reported  this  value  as 
-33.47kcal/mole.  Although  this  last  parameter  may  vary  slightly  with  an 
increase  in  temperature  as  well  as  a  change  in  the  C/metal  ratio,  it  is 
obvious  that  a  large  portion  of  the  total  activation  energy  for  Nb  diffusion  in 
I'lbC^  is  used  for  migration. 

The  similarity  in  values  of  the  '’Nb  diffusion  coefficients  of  the  three 
different  single  crystals  strongly  suggests  that  the  transport  of  the  Nb  in 
the  metal  sublattice  is  uncorrelated  with  the  changing  number  of  C  vacancies. 
Furthermore,  the  structural  arrangement  and  the  octahedral-oct ahedral 
diffusion  mechanism  herein  postulated  for  the  Nb  atoms  in  NbCx  are  equivalent 
to  the  atomic  arrangement  and  transport  process  found  in  elemental  f.c.c. 
metals.  In  this  regard,  Shewman  [17]  has  presented  an  empirical  expression 
for  obtaining  the  activation  energy  Q,  for  self-diffusion  in  f.c.c.  metals. 

This  relationship  is  giver,  by  Q/Tm  36cal/K  where  Tp)  is  the  absolute  melting 
point  of  the  metal.  Using  Tm  =  3873K  for  niobium  carbide,  the  Q  is  found  to 
be  139.4kcal/mole  which  is  in  excellent  agreement  with  the  experimental  valu° 
of  140kcal/mol  obtained  in  this  research. 
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A  comparison  of  the  C  and  Nb-self-diffusion  coefficients  in  all  NbC 
compositions  reveals  that  the  former  species  diffuses  more  rapidly  than 
the  latter  throughout  the  temperature  range  of  this  research.  The  reasons 
for  these  differences  are  primarily  incorporated  in  the  values  of  the 
terms  which  compose  the  activation  energies  for  transport  of  the  two 
atoms.  In  the  case  of  C  diffusion,  in  NbC0.766  the  migration  energy 
necessary  for  this  species  to  move  through  the  large  metal  vacancy  and 
into  one  of  the  numerous  C  vacancies  is  very  small  and  the  value  of  Q  is 
principally  that  of  the  creation  of  the  Mb  vacancy. 

Likewise,  in  MbC0.868  and  NbCo.834,  the  presence  of  the  large  number 
of  structural  C  vacancies  suggests  the  value  of  Q  is  chiefly  derived  from 
the  migration  energy  to  and  from  the  tetrahedral  site  in  the  (0-T-0) 
mechanism.  However,  in  the  Mb  diffusion  mechanism  described  above  the 
activation  energy  term  must  consist  not  only  of  the  energy  of  the  migration 
of  this  species  between  the  C  atoms  but  also  the  energy  of  formation  of  the 
Mb  vacancy. 

B.  High  Temperature  Creep  Studies 
1 .  Equipment  Development 

The  need  for  reliable  creep  testing  of  ceramics  has  necessitated 
the  design  and  fabrication  of  a  new  creep  unit  which  can  apply  compressive 
stresses  which  are  constant  to  within  1%  for  strains  up  to  10%  can  be 
applied  to  the  sample  and  strains  can  be  read  with  an  accuracy  of  5  x  10'7m. 
Furthermore,  loads  as  great  as  440  kg  can  be  applied  to  the  samples  and  the 
furnace  can  be  operated  in  vacuum  or  inert  gas  to  2573K  or  used  with  a 
muffle  tube.  Data  acquisition,  manipulation  and  plotting  is  computer 
controlled.  This  equipment  has  been  described  in  detail  in  ref.  (18),  and 
will  only  be  briefly  outlined  below. 

The  conceptual  plans  of  the  creep  frame  follow  the  fundamental  concepts 
of  Fullman,  Carraker  and  Fisher19  for  the  application  of  constant  compressive 
stress  loading  using  the  Andrade  fulcrum  principle  via  a  lever-arm  mechanism 
from  a  weight  pan  loaded  system.  The  equation 

Li  L2 

!  .  =  _ » 


on 


where  Lj,  L.  and  L3  are  the  distances  between  the  various  knife  edges  shown 
in  Tig.  3,  and  1  is  the  sample  length,  must  be  obeyed  to  maintain  constant 
stress . 


:  I 

! —  t .  ^  .  k~- 


Fig.  3.  Knife  edge  configuration  for  apparatus  used  in  the  present  creep 
studies  of  SiC  and  NbC. 

Two  knife  edges  provide  the  fulcrum  for  the  lever  arm,  another  supports 
the  weight  pan,  and  a  fourth  applies  the  force  to  the  lower  sample  push  rod. 
When  the  load  pan  support  columns  are  vertical,  the  frame  acts  as  a  simple 
fulcrum  to  apply  a  force  on  the  lower  push  rod  which  is  greater'  than  the 
force  on  the  weight  pan  by  a  factor  of  L2/L x .  When  the  sample  deforms,  the 
frame  tilts  backward  and  the  columns  are  no  longer  vertical.  This  movement 
increases  the  force  applied  to  the  push  rod  (and,  therefore  the  sample) 
which  is  proportional  to  the  preset  L3  and  the  amount  of  deformation.  The 
increase  in  force  compensates  for  the  increase  in  area  of  the  sample  as  it 
is  compressed,  thus  maintaining  a  constant  stress.  Shorter  samples  require 
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3  create-'  U  beca-rS  the  area  o*  tho  sa—c’e  beco-e?  create*-  for  t^e  sa-*e 
a-n-’-t  .*  :e*;---a: ;  :r.  t-a-‘  *;•  a  sac'e  :•  .-.•-eate--  ncvr-rt. 

Sar-ple  area  is  not  important  in  r.ak.nn  L,  calculations  because  chanaes  in 
area  for  a  given  amount  of  creep  are  of  the  same  proportion  for  samples 
having  either  a  large  or  a  small  area,  if  both  are  the  same  height. 

The  unit  maintains  constant  stress  to  within  ]"  as  long  as  the  total 
strain  does  not  exceed  10'  and  the  experiment  is  initiated  with  the  weight 
pan  support  columns  in  a  vertical  position.  However,  if  some  idea  of  the 
total  strain  during  the  run  can  be  predicted  and  the  unit  is  set  such  that 
it  passes  through  the  vertical  position  when  approximately  one-half  of  the 
strain  has  occurred,  the  stress  can  be  held  to  within  0 . 5"/  up  to  a  total 
strain  of  10L.  Two  jacuators  which  raise  and  lower  the  main  knife-bearing 
plates  are  used  to  make  this  adjustment  and  are  linked  to  move  simultane¬ 
ously  by  turning  a  single  handle. 

Sample  deformation  is  measured  by  an  LVDT  ultra-precise  gauge  head 
coupled  to  a  signal  conditioner  which  has  a  gauqing  range  of  *5. US  x  10'" m, 
a  repeatability  of  1.016  x  10'7  m,  and  a  linearity  of  '0.2  of  full  scale. 
Since  the  nonlinearity  is  repeatable,  it  is  possible  to  increase  the 
accuracy  of  fhe  creep  measurements  by  performing  a  polynomial  regression 
on  the  calibration  curve  instead  of  a  linear  regression.  Three  curves 
containing  51  points  each  over  the  range  (0  to  5)  x  10'"  m  were  made  for 
the  LVDT  employed  in  the  creep  experiments  using  a  micrometer  calibration 
stand.  The  data  points  were  then  run  through  a  computer  routine  whose 
function  is  to  find  the  polynomial  that  best  fits  the  curve.  However,  the 
sum  of  the  squared  residuals  must  be  at  least  1"’  lower  than  the  next  lower 
order  polynomial  and  must  result  in  a  noticeable  improvement  in  reduction 
of  the  residual  pattern.  For  the  LVDT  system  employed  in  the  present 
research,  a  cubic  polynomial  is  most  satisfactory  and  yields  a  standard 
deviation  of  4.03  x  10"7  m.  This  last  value  is  also  the  smallest  amount  of 
deflection  which  can  be  accurately  measured  with  the  equipment. 

The  previously  mentioned  signal  conditioner  which  powers  the  LVDT  out¬ 
puts  a  dc  voltage  in  the  range  *  1 0  V  which  is  proportional  to  the  LVDT 
deflection.  This  output  is  fed  into  a  dim’tal  dc  voltmeter  which  reads 
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1 0  V  to  a  tenth  of  a  millivolt.  The  voltmeter  has  a  digital  output  and  is 
interfaced  to  a  mini -computer,  which  has  been  programed  to  take  one  voltage 
reading  per  minute.  The  computer  program  is  also  used  to  store  the  data  on 
magnetic  tape  and  to  calculate  true  strain.  The  latter  is  output  to  one- 
channel  of  a  two  channel  flatbed  recorder  via  a  digital  to  analog  converter 
which  is  also  interfaced  to  the  computer.  The  other  channel  of  the  recorder 
is  wired  in  parallel  with  the  dc  voltmeter  to  directly  record  the  voltage 
output  of  the  IVDT. 

After  a  series  of  tests  have  been  completed,  the  computer  calculates 
the  strain  rate,  the  activation  energy,  the  stress  dependence,  and  performs 
any  desired  statistical  analysis.  A  digital  plotter  produces  graphical 
reproductions  of  the  data  in  any  of  a  host  of  possible  formats. 

2.  Creepy  Research 

The  orincipal  investigatcr  has  also  initiated  limited  constant 
stress  compressive  creep.  Studies  have  been  conducted  on  the  same  single 
crystals  of  NbCx  as  those  used  in  the  diffusion  research.  This  research  is 
still  ongoing,  and  a  renewal  proposal  for  its  continuation  is  now  under  review. 
The  results  and  tentative  conclusions  of  this  research,  as  of  this  writing, 
are  briefly  reported  below. 

Samples  0.200  in.  dia.  by  0.300  in.  high  have  been  ultrasonical ly  tre¬ 
panned  from  sections  cut  from  the  lengths  of  the  three  single  crystal  boules. 
The  diffusion  effort  necessitated  the  use  of  only  a  small  amount  of  material 
from  each  boule  and,  therefore,  the  composition  did  not  measurably  change 
from  piece  to  piece  within  a  given  boule.  The  composition  from  one  end  of  the 
boule  to  the  other  does  vary  slightly  (x  =  +  0.02),  however,  as  determined  by 
lattice  parameter  measurements  of  each  of  the  creep  samples.  The  compositions 
noted  below  are  therefore  nominal  compositions  and  are  si ightly.  di fferent 
from  those  noted  in  the  diffusion  discussion. 

To  date,  considerably  time  has  been  consumed  in  learning  to  prepare  pre¬ 
cisely  cut  and  polished  samples  and  in  determining  the  CNb  ratios.  Neverthe¬ 
less,  as  of  this  writing,  ten  creep  experiments  lasting  in  excess  of  1.5 
weeks  each  have  been  conducted  on  samples  previously  annealed  at  2000°C  for 
8  hours  to  remove  residual  stresses.  Such  long  term  experiments  are  required 


because  it  takes  a  very  long  time  (often  in  excess  of  two  days)  for  a  sample 
to  reach  a  steady  state  condition  of  constant  strain  rate  for  a  given  stress 
and  temperature.  The  use  of  on-line  computer  calculations  of  true  strain  and 
strain  rate  and  the  subsequent  storage  of  this  data  allows  one  to  determine 
when  steady  state  conditions  have  been  reached  in  a  given  segment  of  a  run. 
More  information  concerning  the  creep  equipment  and  experimental  procedures 
is  presented  in  the  Plan  of  Research. 

We  have  found  that  all  compositions  undergo  rapid  and  extensive  primary 
creep  even  at  low  stresses  (5-10  x  103  psi)  and  temperatures  (1400  -  1650°C) 
at  the  onset  of  initial  loading.  This  is  repeated  to  a  much  smaller  degree 
for  every  subsequent  increase  in  stress  on  temperature.  As  noted  above, 
steady  state  conditions  do  not  occur,  however,  until  many  hours  (days)  have 
elapsed  from  tne  time  of  application  of  the  new  stress  or  temperature.  This 
could  easily  explain  why  no  previous  investigator  has  established  steady 
state  conditions  in  any  of  the  reported  deformation  or  creep  experiments. 

Most  of  the  creep  research  to  date  has  been  conducted  at  varying  stress 
levels  at  a  given  temperature  in  order  to  determine  the  steady  state  stress 
exponent,  n.  This  value  has  been  determined  to  be  between  4.2  and  4.5  for 
all  three  compositions,  as  shown  in  the  representative  graph.  Fig.  4b,  of 
NbC0.8i  for  T  a  1400°C.  The  corresponding  overall  creep  curve  for  this 
sample  is  shown  in  compressed  form  in  Fig.  4a.  These  values  are  indicative 
of  a  diffusion-controlled  mechanism  of  creep;  however,  a  final  analysis 
awaits  extensive  TEM  on  the  undeformed  and  the  progressively  deformed 
materials.  Equipment  for  thinning  the  many  samples  both  electrochemically 
or  by  ion  milling  is  either  now  established  or  on  order. 

Only  preliminary  research  has  been  conducted  to  date  in  the  determination 
of  the  activation  energy  (Q)  for  steady  state  creep.  However,  .the  extensive 
deformation  of  these  very  pure  materials,  even  at  low  temperatures  and  the 
small  effect  of  increased  temperature  on  the  strain  rate  have  resulted  in 
fairly  low  activation  energy  values.  Although  several  additional  experiments 
are  necessary  to  confirm  these  tentative  results,  it  is  apparent  that  lattice 
self-diffusion  does  not  play  a  role  in  the  steady  state  creep  process. 


1 1 .  Pub  Vicat i ons  Derived  from  th is  Grant 

1.  *J.  D.  Hong,  W.  E.  Griffin  and  R.  F.  Davis,  "Precision  Griridinq 

Device  for  Radioactive  Tracer  Diffusion  Studies  in  Very  Hard 
Ceramic  Materials,"  Rev.  Sci.  Instrum.  49,  83-85  (  1977). 

2.  *6.  B.  Yu  and  R.  F.  Davis,  "Near  Surface  Effect  (NSE)  in  Self- 

Diffusion  in  NbCx,"  Physica  Status  Sol.di  (a)  51,  261  (1979). 

3.  *B.  B.  Yu  and  R.  F.  Davis,  "Thermo-Mechanical  behavior  of  (Jnannealed 

NbC  Single  Crystals,"  J.  Mat.  Sc.  and  Erin.  39,  275  (1979). 

4.  *M.  Hon  and  R.  F.  Davis,  "Diffusion  of  '"C  in  Pol ycrysta  1 1  ine  f.-SiC," 

J.  Materials  Science,  14,  2411  (1979). 

5.  *B.  B.  Yu  and  R.  F.  Davis,  "Self-Diffusion  of  1  C  in  Sinole  Crystals 

of  NbCx,"  J.  Chem.  Phys.  of  Solids,  4C,  997  (1979). 

6.  *C.  H.  Carter,  Jr.,  C.  A.  Stone,  D.  R.  Schaub,  and  R.  F.  Davis, 

High  Temperature  Multi-Atmosphere,  Constant  Stress  Compression 
Creep  Apparatus,"  Rev.  ot  Sci.  Instr. .  5_g,  1352  (1980). 

7.  *B.  F.  Yu  and  R.  F.  Davis,  "Diffusion  of  v,Nb  in  Single  Crystals 

of  NbCx,"  J.  Chem.  Phys.  of  Solids  41,  2088  (1980). 

8.  S.  Chevacharoenkul  and  R.  F.  Davis,  "Hot  Isostatic  Pressing  of 
Niobium  Carbide,"  to  be  submitted  J.  Mat.  Sc.  and  Eng. 


♦Referred  Publications 
I V .  Parti ci oating  Scientific  Personnel 

1.  Robert  F.  Davis,  Principal  Investinator,  Professor  of  Materials 
Engineering,  North  Carolina  State  University,  Raleiqh,  N.  C. 

2.  R.  B.  Yu,  formerly  graduate  research  assistant,  received  Ph.D.  deuree 
May,  1979;  now  employed  as  staff  scientist.  Magnetics  Division  of 
Spang  Industries,  Inc.,  Butler,  PA  16001. 

3.  Calvin  H.  Carter,  Jr.  graduate  research  assistant,  received  M.S. 
degree  December,  1980;  is  currently  working  toward  Ph.D.  degree 
under  the  supervision  of  Robert  F.  Davis. 

4.  S.  Chevacharoenkul,  graduate  research  assistant,  received  M.r.  degree 
May,  1979;  is  currently  working  toward  Ph.D.  deqree  under  the 
supervision  of  Robert  F.  Davis. 


-17- 


V.  Bibliography 

1.  B.  B.  Yu  and  R.  F.  Davis,  Phys.  St  at.  Sol.  (a)  5]_,  261  (1979). 

2.  S.  Sari  an,  J.  Appl.  Phys.  39,  5036  (1968). 

3.  S.  Sarian,  J.  Phys.  Chem.  Solids  33,  1637  (1972). 

4.  S.  Sarian,  in  Recent  Advances  in  Science  and  Technology  of  Materials, 

Vol.  3  (Edited  by  A.  Bishay),  Plenum  Press,  New  York,  p.  247  (1974). 

5.  R.  A.  Andriyevskii ,  Yo.  F.  Khromov,  and  I.  S.  Alekszyeva,  Fiz.  Met. 
Metal! oved  32,  664  (1971). 

6.  N.  N.  Zagriazkin,  Fiz.  Met.  Metalloved.  28,  292  (1969). 

7.  R.  A.  Andriyevskii  and  K.  P.  Gurov.  Fiz.  Met.  Metalloved.  26,  818 

(1968).  — 

8.  VI.  S.  Williams,  Phys.  Rev.  T35,  505  (1964). 

S.  D.  J.  Chadi  and  M.  L.  Cohen,  Phys.  Rev.  B  10,  496  (1974). 

10.  E.  Storms,  B.  Calkin,  and  A.  Yencha,  High  Temp.  Sci .  1^.  430  (1969). 

11.  S.  Sarian,  J.  Appl.  Phys.  40,  3515  (1969). 

12.  t.  Kaufman,  and  E.  V.  Clougherty,  Metallurgy  at  High  Temperatures 
and  High  Pressures  (Edited  by  K.  Gschneider,  N.  Parlee  and  M. 
Hepworth),  p.  322,  Gordon  &  Breech,  New  York  (1964). 

13.  Y.  A.  Chang,  AFML-TR-65-2,  Part  IV,  Vol.  1,  Air  Force  Materials 
Laboratory,  Research  and  Technology  Division,  Air  Force  Systems 
Command,  Wright-Patterson  Air  Force  Base,  Ohio  (1965). 

14.  W.  Schottky  and  C.  Wagner,  Z.  Physik.  Chem.  1_L*  163  (1931). 

15.  L.  Kaufman  and  A.  Sarney  In  Int.Symp.  on  Compounds  of  Interest 
in  Nuclear  Reactor  Technology  (Edited  by  J.  T.  Waber,  P.  Chiotti 
and  W.  N.  Minei),  p.  267.  Met.  Soc.  AIME  (1964). 

16.  E.  Storms,  B.  Calkin,  and  A.  Yencha,  High  Temp.  Sci.  Ij.  430  (1969). 

17.  P.  G.  Shewmon,  Diffusion  in  Solids,  p.  65.  McGraw-Hill,  New  York, 
(1963). 

18.  C.  H.  Carter,  Jr.,  C.  A.  Stone,  R.  F.  Davis,  and  D.  R.  Schaub 
Rev.  of  Sci.  Instrum.  5]_,  1352  (1980). 

19.  R.  1.  Pullman,  R.  P.  Carreker,  Jr,  and  J.  C.  Fisher,  J.  Met.  5^ 

657  (1953). 


